Nematic Quantum Criticality Without Order 
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We consider a two-dimensional interacting Fermi system which displays a nematic 
phase within mean-field theory. The system is analyzed using a non-perturbative 
renormalization-group scheme. We find that order-parameter fluctuations can sup- 
press the nematic order obtained in mean-field theory even at zero temperature. For 
a suitable choice of parameters a quantum critical point surrounded by a disordered 
phase can be realized, giving rise to quantum critical behavior in the absence of an 
ordered regime in the phase diagram. 

PACS numbers: 05.10.Cc, 73.43.Nq, 71.27.+a 

Classical phase transitions are controlled by an energy-entropy balance, and the critical 

behavior of continuous transitions is caused by thermal fluctuations. By contrast, quantum 

l l 

phase transitions occur at zero temperature and are driven by quantum effects [1]. In the 
last decade there has been growing interest in quantum phase transitions in metallic systems, 
and especially in deviations from conventional Fermi liquid behavior due to quantum critical 
fluctuations near a quantum critical point 

as. 

A widely used approach to approximately determine equilibrium phase diagrams of in- 
teracting systems is the mean-field theory (MFT). In numerous cases this approach yields a 
qualitatively correct description of the phase diagram in question. There are however also 
results revealing the inadequacy of this approximation in certain situations, in particular in 



low-dimensional systems. The Mermin- Wagner theorem [4] delivers a well-known, transpar- 
ent and rigorous example of such a result for finite temperatures and spatial dimension d < 2. 
A recent example for the case T = was provided by Greenblatt et al. establishing rig- 
orously the Imry-Ma rounding phenomenon due to disorder for quantum phase transitions. 
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Many other (mostly approximate) results point toward the very important role played by 



quantum fluctuations in d = 2 and T > 



l|-|3) . Also for a number of systems in d = 3 and 



T > (see e.g. Ref. jsj) fluctuations are known to severely influence the phase diagram, for 
example by turning a continuous transition into a first-order one, or vice- versa f|. 

There is considerable interest in two-dimensional itinerant electron systems with inter- 

n n 

actions favoring a <i-wave Pomeranchuk instability [8|, |9[ leading to a so-called electronic 



nematic state 



lOj , where the point group symmetry of the underlying square lattice is spon- 



taneously broken. There is growing evidence that a nematic phase is realized in the ruthenate 
compound Sr 3 Ru207 [U]. Signatures of a nematic instability have also been observed in 
cuprates [12] and other strongly correlated electron materials. 

The nematic transition was found to be generically of first order at zero and sufficiently 



low temperatures for a number of models solved within the mean-field approximation 
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15], while the transition to the normal metallic state becomes continuous only for higher 



temperatures, above a tricritical temperature. In a phase diagram spanned by the chemical 
potential and temperature, the nematic phase resides within a dome-shaped region around 
van Hove filling. Due to the vicinity of the van Hove singularity in the density of states, 
at zero temperature the Landau expansion of the mean-field potential in powers of the 
order parameter exhibits negative interactions at any finite order in 0. Hence, fluctuation 
effects cannot be computed perturbatively. This problem was recently overcome by using 
non-perturbative flow equations derived within the functional renormalization group (RG) 
framework, and it was shown that order-parameter fluctuations may actually suppress the 
first order character of the transition even at T = 0, such that a continuous quantum phase 
transition accompanied by quantum critical fluctuations emerges [r] ]. 

In this letter we show that the impact of quantum fluctuations can be even more drastic. 
We show that the nematic order present within MFT can be completely swept out by the 
order-parameter fluctuations. We also observe that one can tune the system to a scenario 
where a singular quantum critical point is present, while no ordered phase occurs whatsoever. 

We analyse a system of itinerant electrons on a square lattice with a Hamiltonian of the 
form [17] 

H = Y1 Ck ™ k + 2Z /kk'(q) rik(q) n k '(-q) , (1) 

k k,k',q 
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featuring the c?-wave forward scattering interaction 

/kk'(q) = -^(q)dkdk' , (2) 

where d^ = cosk x — cosk y is a form factor with d x 2_ y 2 symmetry. The coupling function 
#(q) > has a maximum at q = and is restricted to small momentum transfers by a 
momentum cutoff. We introduced rik(q) = J2 a c k-q/2 CT c k+q/2 a > wnere c ko- ( c ko-) creates 
(annihilates) an electron with momentum k and spin orientation a, and L denotes the 
number of lattice sites. Note that = nk(Q). We consider a dispersion given by 

ek = — 2t(cos k x + cos k y ) — 4t' cos k x cos k y — 2t"(cos 2k x + cos 2k y ) , (3) 

with nearest, next-nearest, and third-nearest hopping amplitudes t, t', and t", respectively. 

For sufficiently large g = g(0) the interaction drives a <i-wave Pomeranchuk instability 
leading to a nematic state with broken orientation symmetry, which can be described by the 
order parameter 

*=££><n k >. ( 4 ) 

k 

Introducing a fluctuating order parameter field via a Hubbard- Stratonovich transformation, 
the fermionic interaction term in Eq. ([1]) is decoupled, and the fermionic degrees of freedom 
can be integrated out. Retaining only the leading momentum and frequency dependencies of 
;he two-point function and neglecting such dependencies in the higher-order vertex functions 



181 ] . one obtains the following action: 



where q ,a; n is the momentum representation of the order parameter field 0, and u n = 2irnT 
with integer n denotes the (bosonic) Matsubara frequencies. The momenta and frequencies 
contributing to the action S[<p] are restricted by an ultraviolet cutoff A to the region y4 -^ + 
Z q 2 < Z Aq. In the fermionic representation, A is related to the maximal momentum 
transfer allowed by the interaction in Eq. (JTJ. The values of the prefactors A and Z are 
related to the <i-wave particle-hole bubble in presence of a (mean-field) order parameter 6$ 



and Z depends also on the specific form of the function g(q), as discussed in Refs. [la. Il9|. 
We do not treat a possible dependence of A and Z on the direction of q, such that the 
values we insert should be viewed as an angular average. The effective potential U[<j>] is 
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given by 

U[(f>]= j dr I d 2 rU((f)(r,T)) , (0) 



T 



where 

U{4>) = f g - 2T J In (1 + e -^-^' T ) . (7) 

An action of the form (jSJ) was already considered in Ref. 16J. A polynomial expansion of 
U(4>) yields negative coefficients at any order in cf) and is therefore inadequate here. Our 
calculations are carried out within the one-particle irreducible version of the functional 
RG [20]. The analysis amounts to solving the flow equations for U(<p) and Z derived in 
Ref. 16j. Here we only remark that we employ a truncation of the functional RG based 
on the approximation that the effective action retains the local structure given by Eq. (JBj). 
We do not impose any ansatz for U(4>) and therefore keep track of the renormalization of 
infinitely many couplings. The flow of the coefficient Z captures the anomalous dimension 
of the order parameter field at criticality. The truncation level is analogous to the next-to- 
leading order in the derivative expansion, which was developed in the context of classical 
criticality |2l|,|22|. 

We now present our results, which are obtained from the solution of the RG flow equations 
for the following set of parameters: t = 1, t' = —1/6, t" = 1/5, g = 1/2, A = 1, Z = 10, 
and for different values of A . For this choice of hoppings the density of states is strongly 
asymmetric and the chemical potential corresponding to the maximal T c is separated from 
the van Hove energy (/i V H = —1.444). The relatively large value of Z corresponds to a 
pronounced peak of g(q) at q = 0, which favors the nematic instability compared to insta- 
bilities at a finite q. Fluctuations reduce the magnitude of the order parameter compared 
to the mean-field value, as expected, and the nematic region in the phase diagram shrinks 
(see Fig. 1). The suppression is naturally stronger for larger cutoffs A . For sufficiently 
large A the nematic phase is located to the right of the van Hove filling. When A exceeds 
a critical value Aq, no ordering occurs whatsoever. By continuity, for Ao = Aq the phase 
diagram exhibits an isolated quantum critical point at /i = //* and T = 0, while the system 
is disordered in the entire (fi,T) plane. In the vicinity of the point (/x = /iq C ,T = 0) the 
system displays soft fluctuations characterized by a large correlation length. For the present 
choice of parameters we find /x* c ~ —1.404 and Aq w e -°- 2574 . Note that for smaller values 
of g the order would be destroyed already at smaller Aq. 
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FIG. 1: (color online). Critical temperature for the nematic transition as a function of the chemical 
potential /x. Left panel (a) Critical temperature in mean-field theory (MFT) and in the presence 
of fluctuations with four different choices of the ultraviolet (momentum transfer) cutoff, Ao = 
e -1 , e -0 - 4 , e -0 - 28 , e -0 - 26 . Right panel (b) Enlargement of the results for Ao = e -0 - 28 , e~ a26 . 

We have also computed the behavior of the correlation length £ upon approaching the 
quantum critical points at the edges of the nematic dome for Ao < Aq, and for the isolated 
quantum critical point in the special case A = Aq. The correlation length in the symmetric 
phase is determined by the curvature of the potential U"(0) and the renormalization constant 
Z at the end of the flow via £ 2 = Z/U"(0). For A < Aq the results are consistent with 
the predictions of Hertz-Millis theory [23] for a two-dimensional system with a dynamical 
exponent z — 3, namely £ 2 ~ |TlogT| _1 at fixed \i = /i qc and £ 2 ~ — /iqd -1 at T = 0. 

The scenario occurring for A = Aq requires fine-tuning. However, features resembling 
quantum-critical behavior persist also for Ao > Aq. In Fig. 2 we plot the correlation length 
£ as a function of fi at T = for a sequence of values of A approaching Aq from above. 
The peak formed around \i = /i* becomes rapidly pronounced as A — > Aq. Remarkably, 
although the quantum critical point occurring for Ao = Aq is Gaussian, the exponent v 
describing the divergence of £ at // — > /jl* takes the value 1 instead of the usual value ~. 
This is easily understood by noticing that for A = Aq the quadratic mass term in the 
Landau expansion must not take negative values in the absence of an ordered phase and 
therefore cannot be proportional to /i — fi* as in the usual cases. The natural /^-dependence 
of the quadratic term is therefore quadratic in \l — \f near //* which leads immediately to 
v — 1. This peculiarity does not influence the behavior when the singularity is cut off by 



temperature. In fact for fixed \i = fi* we recover the usual 23J behavior £ ~ |TlogT| 1 / 2 , 
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FIG. 2: (color online) Top (a) Correlation length £ as a function of jx at T = for various values 
of Aq > Aq. Bottom (b) Double logarithmic plot of £ versus /i — /x* c and fit to the exponent v = 1. 



as shown in Fig. 3. For A > Aq no singular point in the phase diagram occurs, but the 
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FIG. 3: (color online) Temperature dependence of the correlation length £ for the special case /U = 
fj,t c and Ao = Aq. Inset: Double logarithmic plot exhibiting the expected behavior £ ~ l^log T| -1 / 2 
for small T. 

correlation length takes sizable values in the vicinity of the point (fi = /iq C ,T = 0) as long 
as Aq is not significantly larger than Aq. 
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One may be surprised to find the nematic order fully eliminated by fluctuations even at 
van Hove filling. In MFT, the diverging density of states at van Hove filling leads to an 
instability of the symmetric phase at any (arbitrarily small) coupling g 15]. One would 
expect this to remain true also beyond MFT, provided the system remains a Fermi liquid. 
Inversely, the absence of nematic order at van Hove filling despite the presence of a <i-wave 
attraction g would imply the absence of fermionic quasi-particles, that is, non-Fermi liquid 
behavior. Unlike the situation at the quantum critical point, that non- Fermi liquid behavior 
is not related to critical fluctuations. Working with an effective order parameter action we 
do not have direct access to fermionic properties. However, the breakdown of Fermi liquid 

n 

theory at van Hove points was already diagnosed many years ago by Dzyaloshinskii [24|. In 
the light of that work the complete suppression of order by fluctuations even at van Hove 
filling is a consistent possibility. 

In summary, we have analyzed a two-dimensional interacting electron system which ex- 
hibits a nematic phase in mean-field theory. We have shown that fluctuations can extinguish 
the nematic order completely from the phase diagram spanned by the chemical potential 
and temperature. For a special choice of parameters a quantum critical point can be re- 
alized in the absence of order. Despite being realized as a consequence of strong quantum 
fluctuations, the quantum critical point itself is Gaussian. The temperature dependence of 
the correlation length in the quantum critical regime above this isolated quantum critical 
point follows the conventional behavior known for Gaussian quantum criticality. However, 
at zero temperature the correlation length diverges with an exponent v = 1 instead of v = | 
upon approaching the critical chemical potential. While the present calculation was carried 
out for a specific interaction triggering a nematic phase, an analysis of fluctuation effects 
in other two-dimensional systems described by an action of the form (jSJ) can be performed 
along the same lines. 
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